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The  actual  choice  of  a  powered  flight  trajecto^^  foa  the  farth- 

comiiig  lOY  wtll^dependSupon  details  of  booster  design  and  upon  the 

particular  orbit  selected,  the  latter  being  determined  by  instrumentation  ar.d 

data  requirements.  Independent  of  these  considerations,  however,  it  is  of 

interest  to  ask  what  type  of  trajectory  will  be  optimal  from  the  standpoint 

of  missile  efficiency,  ^he  resets  of  such  an  investigation  are  reported  here. 

U  f  0^ 


To  begin  with  we  must  select  some  criterion  of  excellence  for 
comparing  trajectories.  Whereas  range  is  an  obvious  choice  for  surface  to 
surface  missiles,  it  is  not  clear  what  quantity  will  play  a  similar  role  for 
satellites.  One  important  property  of  a  satellite  is  its  useful  orbiting  life¬ 
time.  Too  low  an  orbit  altitude  will  result  in  a  rapid  slowing  down  due  to 
atmospheric  drag  and  an  undesirably  short  lifetime.  Thus,  if  the  perfor¬ 
mance  of  the  booster  vehicle  is  marginal,  orbit  altitude  will  be  of  crucial 
importance  and  can  serve  as  the  desired  criterion.  This  may  not  be  a 
sensible  choice  if  the  booster  is  sufficiently  powerful,  for  at  altitudes  above 
about  300  miles  the  orbiting  life  time  is  measured  in  years  and  it  is  likely 
that  the  power  supply  life  rather  than  the  drag  induced  decay  would  be  the 
determining  factor.  On  the  other  hand,  in  such  a  case  considerations  of 
missile  efficiency  are  anyhow  of  littl.»  importance. 

The  orbit  "altitude"  is  a  well  defined  quantity  only  for  circular 
orbits.  In  the  case  ox  elliptical  orbits  it  is  convenient  to  work  with  the 
distance  of  closest  approach  (d.c.a.),  i.e.,  the  altitude  of  the  satellite 


(1)  H.  E.  Newell,  Scientific  Uses  of  an  Attifif  ial  Earth  Satellite,  Jet  Propulsion 
25  717  (1955) 
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at  perigee.  This  is  approximately  correlated  with  the  drag  induced  decay, 
since  the  major  part  of  the  slowing  down  will  occur  near  perigee  The 
correlation  is  not  exact,  for  of  two  orbits  with  equal  d.c.a.  the  one  with 
the  greater  eccentricity  will  have  the  longer  life.  However,  the  analysis 
described  below  would  not  be  changed  essentially  if  the  d.c.a.  were  re¬ 
placed  by,  say,  the  mean  of  the  altitudes  at  perigee  and  apogee  or,  in 
general,  by  any  function  of  the  altitude  and  vector  velocity  at  the  end  of 
powered  flight .  For  definiteness,  we  shall  adopt  the  d.c.a.  as  our 
criterion. 

Analytical  results  on  trajectories  can  only  be  obtained  by 
neglecting  certain  terms  in  the  equations  of  motion,  for  if  we  include 
atmospheric  effects,  variations  in  the  gravitational  force  field,  etc.,  then 
solutions  of  the  differential  equations  can  only  be  obtained  by  numerical 
integration.  It  is  rather  gratifying  to  find  that  with  a  few  ,  not  unreasonable 
approximations  we  are  led  to  a  very  simple  result  the  tangent  of  the  thrust 
attitude  angle  should  be  a  linear  function  of  time.  In  the  following  section 


the  derivation  of  this  result  and  the  method  of  determining  the  coefficients 
in  the  linear  relation  will  be  presented.  We  shall  see  that  while  this  result 


y  is  an  approximate  one,  it  can  be  of  practical  value  when  correctly  applied 
-OiM'  approach  wtil  ba  assumeithat  the  characteristics  of  the 


booster,  i.e.,  its  thrust  and  mass  as  functions  of  time,  are  specified.  The 


problem  is  then  to  determine  what  powered  flight  trajectory  will  put  the 


satellite  into  an  orbit  of  maximum  altitude, 
are  possible,  e.g.:  for  given  payload  (i.e. 


,  orln 


any  alternative  formulations 


iting  weight)  and  orbit,  find 


the  trajectory  which  minimizes  the  takeoff  weight.  Problems  of  this  sort. 
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however, ,  require  fairly  detailed  Information  about  the  booster  --  staging 
configuration,  atructural  efficiencies,  tankage  factors,  etc.  »bo  that 
solutions  of  general  validity  can  scarcely  be  obtained.  These  difficulties 
are  avoided  when,  as  here,  the  properties  of  the  missile  are  regarded 
as  given. 
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II.  TRAJECTORY  ANALYSIS 

Consider  the  problem  of  fir.duig  the  powered  flight  trajectory 
which  will  result  in  maximum  satellite  orbit  altitude  for  a  missile  whose 
thrust  magnitude  and  mass  are  specified,  albeit  arbitrary'^^  functions 
of  time,  F(t)  and  M(t).  We  shall  treat  the  missile  as  a  point  mass,  the 
earth  as  a  non-rotating  sphere,  and  to  begin  with  shall  neglect  three 
things:  (1)  Aerodynamic  effects  -  drag,  lift.  etc. 

(2)  The  dependence  of  thrust  on  altitude 

(3)  The  variation  of  the  magnitude  and  direction  of  gravity 
during  powered  flight 

If  the  orbit  is  to  be  circular  (elliptical  orbits  are  discussed 
later)  then  at  booster  burnout  the  vertical  component  of  velocity  must  be 
zero  and  the  tangential  component  must  have  such  a  value  that  the  centripetal 
and  gravitational  forces  just  balance, 


1/2 


iiat 


=  (gR^I  (I  t  y/R^) 


71/2 


u 


(1) 


where  (u  .  u  )  is  the  booster  velocity  at  burnout,  R  is  the  earth's  radius, 

X  y  e 

y  is  the  burnout  altitude  and  g  is  the  acceleration  of  gravity  at  the  earth's 
surface.  For  planar  trajectories,  the  path  is  completely  determined  by 
specification  of  the  thrust  attitude  angle  (say  with  the  horizontal),  Ij./  (t). 

We  then  ask  'vhat  function  'y  (t)  maximizes  the  altitude  at  burnout  subject 
to  the  conditions  (1).  For  simplicity,  we  shall  first  assume  to  have  a 

given  value,  ignoring  its  dependence  on  altitude.  Later  we  shall  sea  that  thus 
dependenr”  can  easily  be  taken  into  accoant. 
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If  6y  it  the  change  in  burnout  altitude*  dup  to  a  change 
(t)  +  6  with  similar  definitions  for  6u^  and  6Uy, 

then  a  necessary  condition  for  y  to  be  a  maximum  subject  to  ( 1)  is 


6y  -  X6u  -  ^6u 
X  y 


(2) 


where  the  Lagrange  multipliers  k,  p  are  unknown  constants.  If  N(t) 
=  F(t)/M(t)  is  the  instantaneous  thrust  to  mass  ratio,  we  have 


u 


pT 

I  dt  N  cos 

<Jo 


n 


dt 


N  sin  ij/  -  gT 


o 


f 


dt  (T  -  t)  N  sin  y  -  1/2  gX' 


(3) 


where  T  is  the  total  burning  time.  If  6  {t),  then 


u  — ^  u 

X  X 


dt  N  sin  y  6  jC(t) 


Jo 


so 


6u 


n  T 


O 


dt  N  cos  y6  y 


jp,p, 


(4) 
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cos 


H) 


T 

^  dt  (T  -  t  )  N  co8(^ 


Upon  substituting  these  expressions  into  (2),  we  see  that  (2)  will  be 
satisfied  for  arbitrary  6'L^(t)  if  and  only  if 

N  ^(T  -  t)  cosl|*+  X.  sinl^  -  p  co8\|.^|  =  0 
for  all  t.  Thus,  when  N  /  0  we  require 


tan 


*  T)  +  f| 


Since  \  and  p  are  anythow  unknown  constants,  we  may  simply  state  our 
result  in  the  form 


tan 


a  -  bf 


(51 


where  a  and  b  are  constants  which  must  be  chosen  to  satisfy  (1). 

Although  the  result  (5)  is  valid  lor  arbilr^_y^^^  F(t)  and  M(t),  the 

explicit  values  of  a  and  b  will,  of  course,  depend  upon  these  functions.  In  the 

particular  case  where  F(t)  and  M(t)  are  (piecewise)  constant  functions,  the 

expressions  (3)  for  u  and  u  can  be  integrated  in  terms  of  elementary  functions 

X  y 

The  increments  in  u  ,  u  and  y  due  to  one  stage,  of  mass  ratio  r,  exhaust 

X  y 


velocity  c  and  burning  time  tj,  when  tan  is  a  linear 


function  of  time  are 


(2)  The  term  "arbitrary  function"  is  to  be  understood  in  a  physical  rather  than  a 
mathematical  sense,  i.e.,  as  a  function  which  is  arbitrary  except  for  such  re¬ 
strictions  of  continuity,  differentiability,  positive  definiteness,  etc.,  as  may  be 
necessary  to  exclude  situations  which  are  clearly  not  sensible  in  view  of  the 
physical  significance  of  the  function. 
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Au  = 

X 

c  • 

c  o  9  , 

^  Inr  - 

>„  , 

1  f  c  08 

1 

1  t  C  08 

'Vi 

■  ^) 

Au  = 

y 

Au 

X 

tan  + 

cln  1 

1  — 

•‘"Vo 

CO.  y, 

I- 

sin 

cos 

Ay  = 

‘‘  1 

c  cos 

v. 

•  cos 

Vo 

Au 

y 

\ 

sin 

) 

o 

V 

T  -  1 

(6) 


where  is  the  initial  value  of  is  the  value  at  time  t^  aud  is 

the  value  at  time  t  =/m  \  =  t^r  ,  M  and  M  being  the  initial  mass  and 

Att)  “ 

the  mass  flow  rate  for  the  stage  Bv  summing  the  contributions  of  each  stage 

and  adding  (to  u^  and  y)  the  gravitational  terms  (-gT  and  -1/Z  gT^  ,  where  T 

is  the  total  burning  time),  we  obtain  expressions  for  u  ,  u  and  y.  Setting 

X  y 

u  =  u  .  and  u  =0,  we  would  use  two  of  these  relations  to  find  a  and  b 

X  sat  y 

and  then  substitute  the  resulting  vaiues  of  a  and  b  into  the  expression  for  y  to 
find  the  orbit  altitude. 

Although  we  have  demonstrated  »  nly  that  (5)  is  a  stationary  solution, 
we  note  from  (3)  and  (4)  that 

T 


6^y  = 


dt  ( T  -  t  )  N  sin  (6  t|>') 


will  be  negative  for  arbitrary  6  ’^-’provided  siny  ^  0.  Thus,  we  are  guaranteed 
that  y  is  actually  a  maximum  provided  that  a  and  b  satisfy 
a  -  bT  0 


(3)  In  terms  of  the  a  and  b  of  (5),  tan  4^^-  ^  ^  tan4^j  =  a-btj^  tan'^^ 


a  - 


(bt  j  r 

V~7-t 


i 
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While  the  aclu.il  jmih  ''S';  *l«'terming  a  and  b  is  straight  forward, 

it  is  somewhat  lengthly  in  |ir-ii  1 1«  i  .  .-.imc  it  involves  the  solution  of  two 

simultaneous  tran8cendenf«il  oqu.itions.  Also,  our  derivation  has  neglected  the 

dependence  trf  u  upon  altitudn  IJolh  of  these  difficulties  can  be  resolved  if 

we  adopt  a  somewhat  mori-  m  ■'  ii  imo.uIi 

We  drop  the  previous  restriction  lo  circular  orbits  and  simply 

assume  that  at  burnout  we  h.ive  i  eiiam  values  of  u  ,  u  and  y.  These  uniquely 

X  y 

determine  the  Kepler  ellip:-'*  whieli  Mu:  salellile  follows  and,  in  particular,  the 
distance  of  closest  appro.u  fi  {d.i  a  )  of  this  oil>it.  If  H  's  the  d.c.a.  and 


D  =  +  H 

e 


then  in  terms  of  energy,  E,  €ind  angul.ir  rnomenf  ri,  Li,  we  have 


where 


D  =  -k  (k^  t 
ZE 

E  =  u^  -  k(R^.  +  y)’  ‘ 


(7) 


L  =  u  (R  +  vl 
X  e 


k  =  gR, 


(8) 


If  we  simply  choose  L^(t)  1o  maximjre  D.  then  we  require 
6D  =  D  6y  +  D  6u^  +  6u  =  0 

^  X  V 


(9) 


where  D^,  =  3D/3y.  etc.  Using  (4)  we  again  find 

tan  ^  =  a  •  bt 
where  now  a  and  b  are  given  by 

a  =  (TD  ^  )/D^  b  =  D  /D^ 

v  X  X 


(10) 


I 


t 


I. 


I 

'I 


It 
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Thie  provides  a  systematic,  iterative  procedure  for  calculating 

a  and  b:  guessu  ,  u  ,  y;  calculate  D  ,  D  ,  D  from  (7)  and  (8);  compute 
xy  uuv 

X  y  ^ 

a  and  b  from  (10);  find  u  ,  u  ,  y  from  (6);  etc.  This  will  yield  the  maximum 

X  y 

y  possible  with  this  missile,  subject  to  the  assumptions  (1),  (2),  and  (3). 

Of  these,  the  most  serious  is  the  first,  since  a  program  like  (5) 
would  result  in  intolerably  large  angles  of  attack.  Consequently,  it  seems 
advisable  to  use  a  "gravity  turn"  (thrust  parallel  to  velocity)  until  the  missile 
is  high  enough  to  justify  the  neglect  of  aerodynamic  forces,  at  which  time  we 

(4) 

can  switch  over  to  the  program  (5).  Using  the  fact  '  '  that  the  velocity  magnitude 
at  the  end  of  a  gravity  turn  la  a  far  less  sensitive  function  of  the  initial  condi¬ 
tions  (kick  angle,  etc.  )  than  is  the  velocity  •tti’ude,  it  is  easy  to  show  that 
the  gravity  turn  should  be  chosen  so  that  a*  the  time  of  transition  to  the 
regime  (5),\|^  has  no  discontinuity. 

This  prog  ram  also  takes  care  of  assumption  (2),  since  we  follow  (5) 
only  at  high  altitudes  where  the  thrust  is  independent  of  alfitu'i'-  Finally,  the 
fact  that  the  powered  flight  takes  place  in  a  non -constant  gravitational  force 
field  would  modify  (23)  somewhat.  However,  this  effect  will  be  small  pro¬ 
vided  the  distance  covered  during  powered  flight  is  not  large  compared  to  the 
earth's  radius. 


(4)  G.  Culler  and  Burton  D.  Fried,  'Universal  Gravity  Turn  Curves'', 
to  be  published 
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IV.  CONCLUSIONS  AND  DISCUSSION  OF  RESULTS 

We  have  shown  that  to  obtain  a  satellite  orbit  of  maximum  altitude 
with  a  given  missile,  the  thrust  attitude  angle,  1^,  should  vary  with  the  time 


according  to 


a  -  bt 


during  that  portion  of  the  powered  flight  where  aerodynamic  forces  and 
variations  in  gravity  are  unimportant.  This  result  is  valid  for  arbitrary 
time  dependence  of  thrust  and  mass,  and  hence  for  arbitrary  staging  configu¬ 
rations.  Moreover,  if  we  replace  orbit  altitude  by  any  other  function  of  the 
burnout  altitude  and  velocity  vector,  e  g.  the  altitude  at  perigee,  and  look 
for  a  stationary  solution,  we  again  obtain  a  result  of  the  form  (5). 

The  evaluation  of  a  and  b  requires  that  we  calculate  the  burnout 
velocity  with  given  by  (5)  For  the  particular  casi?  where  thrust  and 

mass  flow  rate  are  constant,  this  can  be  carried  through  analytically.  The 
results  -  eqs.  (6)  -  are  a  generalization  of  the  familiar  clnr  formulas 

which  are  obtained  when  is  constant 

A  steering  program  like  (5)  is  physically  reasonable.  To  obtain 
high  altitudes,  we  would  expect  to  start  firing  with  large  ,  so  that  the 
resultant  vertical  velocity  could,  in  the  course  of  time,  produce  altitude. 
Later  in  the  flight  we  would  want  smaller  in  order  to  build  up  the 

required  horizontal  velocity  and  also  to  allow  gravity  to  cancel  any  accumu¬ 
lated  vertical  velocity.  The  optimum  way  of  combining  these  considerations 
is  prescribed  by  (5).  Since  the  general  form  of  this  function  is  qualitatively 
similar  to  that  associated  with  a  gravity  turn,  we  can  expect  that  the  penalty 
of  using  the  latter  during  the  atmospheric  portion  of  the  flight  will  not  be 
large. 
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